The role of carbon for superconductivity in MgC x Ni 3 from specific heat 
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The influence of carbon deficiency on superconductivity of MgCNi 3 is investigated by specific heat 
measurements in the normal and superconducting state. In order to perform a detailed analysis of the 
normal state specific heat, a computer code is developed which allows for an instantaneous estimate 
of the main features of the lattice dynamics. By analyzing the evolution of the lattice vibrations 
within the series and simultaneously considering the visible mass enhancement, the loss in the 
electron-phonon coupling can be attributed to significant changes of the prominent Ni vibrations. 
The present data well supports the recently established picture of strong electron-phonon coupling 
and ferromagnetic spin fluctuations in this compound. 

PACS numbers: 74.25.Bt, 74.25. Kc, 75.40.Cx 



The discovery of superconductivity in MgCNi 3 has 
caused much attention^ since the large Ni content sug- 
gests a magnetic state rather than superconductivity. In- 
deed, up to now a lot of experimental and theoretical 
publications point to a ferromagnetic instability at tem- 
peratures near and below the superconducting transition 
temperature T c ss 7 KiSiMaiii So far there is much discus- 
sion on the carbon content in MgCNi 3 , since only sam- 
ples with carbon excess of ss 50 % show large T c values, 
whereas for the stoichiometric composition a strongly re- 
duced T c is found?i& which may be triggered by enhanced 
pair-breaking due to increasing spin fluctuations. 9 How- 
ever specific heat measurements indicate that this is not 
the case but instead the electron-phonon coupling is sig- 
nificantly reduced^ accompanied by a considerable hard- 
ening of low-energy vibrations^ Band structure calcula- 
tions predict a constant 9 or decreasing^ electron density 
of states at the Fermi level (EDOS), leaving the interest- 
ing prospect of significantly changing lattice dynamics. 
Recent measurements of the carbon isotope effect sup- 
port this pictured 

Polycrystallinc samples of MgCNi 3 have been prepared 
by solid-state reaction. The high volatility of Mg was bal- 
anced in the usual way by Mg excess. Three samples of 
nominal composition Mg^^ 2 C Xn Ni3 with x n = 0.75, 0.85 
and 1.00 (n=nominal) have been prepared. For compar- 
ison the previously reported sample with x n = 1.60 and 
T c — 6.8 K is chosen^ Details on the preparation proce- 
dure have been published elsewherei The obtained sam- 
ples were characterized by x-ray diffractometry. Due to 
the relatively low carbon content no additional graphite 
was found in the samples, instead a small fraction of 
MgNi 2 forms with decreasing carbon content. However 
the fraction never exceeds 6 Vol.-% and has minor influ- 
ence on the following specific heat analysis. The positions 
of the reflexes in the x-ray diffractograms are shifted to 
higher angles for decreasing carbon content. The left 
panel of Fig. ^ shows a magnification of the region 
20 « 84 deg, where the influence of the carbon deficiency 
is visible. For a full diffractogram of the x n = 1.60 sample 
see Ref. |(J The lattice constants of the samples were de- 
termined using the Rietveld code fullprof. 12 From that 



the effective carbon content was estimated according to 
Ref. 0- The dependence is plotted in the right panel 
of Fig. ^ together with the estimated superconducting 
transition temperatures^ It is obvious that the present 
samples with x n < 1.00 are multi-phase samples. How- 
ever this is of minor importance for the following analy- 
sis. The specific heat was measured for T = 2 — 300 K 
and magnetic fields up to fioH = 12 T using a Quan- 
tum Design Physical Property Measurement System. For 
the sample with x n = 1.00, which shows a considerably 
broadened superconducting transition (due to its multi- 
phase nature) the measurement was extended down to 
T = 0.3 K. 



I. COMPUTATIONAL DETAILS 

From the mathematical point of view, the extraction of 
the phonon density of states (PDOS) from specific heat 
measurements is ill-posed. This means, that a number of 
very different PDOS can lead to very similar specific heat 
curves. From the experimental point of view one has to 
extract the phonon specific heat from the electron and 
other non-phononic background exactly. In reality tem- 
perature dependent effects (e.g. softening) then prevent 
an inversion approach in most cases. From this start- 
ing point we developed a computer code based on simple 
models, namely a Debye and an Einstein model, in order 
to get a rough estimate on the lattice dynamics from spe- 
cific heat measurements instead of a mathematical exact 
inversion. Specifications of the code, which is fully based 
on ROOTjH will be published elsewhere. Here we want 
to sketch only the main ideas and illustrate the capabil- 
ities of the code. The specific heat of a compound with 
n atoms per unit cell has 3n — 3 optical and 3 acoustic 
vibrations which can be approximated by a linear combi- 
nation of Einstein and Debye models. The contribution 
of the Einstein phonons to the specific heat is given by: 
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Figure 1: Structural influence of nominal carbon content on 
MgC^Nia. Left panel: x-ray diffractogram around 20 « 
84 deg. Right panel: relation between lattice constant a, ef- 
fective carbon content x e and T c according to Ref. for the 
investigated samples. 



Figure 2: Energy resolution AOoi/ODi and A0Ei/6>Ei of the 
code in a background of 6 (5) random Einstein and 1 (2) ran- 
dom Debye terms. Left panel: energy resolution for one Ein- 
stein vibration. Right panel: energy resolution for one Debye 
vibration. 



with number of optical modes N and Einstein tempera- 
tures Oei- The Debye contribution reads: 
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with Debye temperatures Qui- The form of the phonon 
density of states is thereby simplified as: 
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with step function 9 (xq — x) and the characteristic tem- 
peratures in meV, indicated by and wgi- For the 
case of MgCNi 3 one is left with 3 Debye and 12 Einstein 
terms, adding up to 15 parameters. To avoid the need of 
restrictions on the parameters we developed a fast con- 
verging algorithm. It makes use of a library containing 
integrated values of a single Debye and Einstein model 
for AT = 1 K temperature intervals between T = 2 K 
and 300 K and the overall integral of the two models in 
the interval T = 2 — 300 K. The algorithm integrates the 
measured specific heat in a number of temperature inter- 
vals and once over the measurement range T = 2 — 300 K. 
The use of the overall integral for the fitting procedure 
accelerates the procedure considerably and most impor- 
tantly allows the omission of parameter-restrictions. The 
code was statistically tested by randomly distributing 6 
Einstein and 2 Debye terms. The fitting result was than 
compared to the starting conditions. Fig. [21shows the re- 
sulting energy resolution gained from 20000 analyzed spe- 
cific heat model-curves. This test shows that the energy 
resolution is weak for low-energy optical and high-energy 
acoustic vibrations, which are both rare cases. However, 
there is still potential in the algorithm to improve the en- 
ergy resolution. For example by using a set of weighting 
factors for different temperature regions, which we are 
currently working on. 



II. ANALYSIS AND DISCUSSION 

The specific heat of the sample with x n = 1.00 is shown 
in Fig. |31 (the samples with x n — 0.85 and 0.75 show very 
similar results and are omitted for clarity). The specific 
heat of MgCNi 3 is given by a lattice part entice (T), de- 
termined by the above mentioned code and an electron 
part 7nT with the Sommerfeld parameter: 



7n = [1 + A ph + A sf (0)] 70, 



(1) 



with free electron parameter 70 = 11 mJ/ (molK 2 ). The 
mass enhancement due to the electron-paramagnon cou- 
pling can be estimated from 
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where tp (iz) is the digamma function and z = 
lj j (27rfcBT). The electron-paramagnon spectral density 
is given by 

a 2 F s f (lj) — auj9 (lj s { — lj) + 9 (lj — lj s {) ■ 

The paramagnon energy, which was estimated in the case 
of x n = 1.60 as lj s { w 2.15 meV?& was also used for the 
present samples. The black line in Fig. 131 is the result of 
the fitting procedure. The differences between the data 
and fit results are shown in the inset for x n = 0.75, 0.85 
and 1.00. They are very small, in particular for x n = 0.85 
and 1.00, indicating the success of the fitting procedure. 

Fig. 01 shows the low temperature region of the spe- 
cific heat. The reduction of 7n and the increasing flat- 
tening of the curves with decreasing car bon content was 
already discovered in a previous work bv lShan et alfe For 
x n = 1.00 a strongly broadened superconducting transi- 
tion is visible at T c « 2.5 K, in accord with the mentioned 
multi-phase nature. The deviation from the low temper- 
ature Debye approximation (linear behavior) due to spin 
fluctuations is clearly seen for all samples (Fig. |3J|. 
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Figure 3: Specific heat of MgCj 00 Ni3. Solid line: specific heat 
fit (see text for details). Inset: relative difference between fit 
and specific heat for x n = 1.00, 0.85 and 0.75 (from top to 
bottom). 



Figure 4: Specific heat of MgC^Nis at low temperatures. 
Solid lines: fits to the data (see text for details). The anomaly 
at T m 2.5 K for sample x n = 1.00 indicates the supercon- 
ducting transition. Inset: entropy change of the electronic 
specific heat within the superconducting state for x n = 1.00. 



The model description of the normal state is given as 
solid lines. For x n = 1.00 the conservation of entropy of 
the superconducting state, which is expected for a phase 
transition of second order, was used as an additional re- 
quirement for a successful fit (inset of Fig. 7n was 
varied with a step-size of 0.05 mJ/molK 2 until the rela- 
tive difference shown in the inset of Fig. 0|was minimized. 
Close to the final value of 7n, the derived sets of vibra- 
tional parameters did not appear qualitatively different 
from each other. The results of the fitting procedure 
are summarized in Tab. [I] The electron-paramagnon 
coupling constant A s f (0) only slightly decreases by low- 
ering the carbon content, despite a predicted increase 
of spin fluctuations & The electron-phonon coupling con- 
stant shows a strong decrease from A p h = 1.84 (x n = 
1.60) down to Aph = 1.06 (x„ = 0.75), in agreement with 
results obtained bv lShan et alfe 

Fig. [5] shows the resulting phonon density of states, 
F (lu) for x n — 0.75, 0.85 and 1.00 in comparison with 
previous results for x n = 1.60.- There are two main 
changes visible within the series. First there is a consider- 
able change of the high-energy mode, which is expected, 
since it is dominated by carbon. Second all low-energy 
modes are slightly shifted to higher energies, accompa- 
nied by the shrinking of the unit cell. The extent of this 
behavior can be quantified by calculating the character- 
istic phonon frequency: 

[ 2 f°° a 2 (ui)F(u) 1 
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and a 2 (w) as the electron-phonon interaction function. 
In a previous work we suggested to use an approach 
of the form a 2 (u>) — S/y/uj, with scaling parameter 
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Figure 5: Phonon density of states of MgC^Nis. Dashed lines: 
result for x n = 1.60£ 



5 to approximate the electron-phonon interaction func- 
tion. Within this approach the low-energy phonons are 
more strongly weighted. Calculating u)\ n for the present 
samples, an increase from u>\ n = 154 K (x n = 1.00) to 
162 K (x n = 0.85) and 164 K (x n = 0.75), compared to 
LU\ n = 143 K for x n — 1.60 is founds Even by leaving 
out the carbon-dominated high-energy mode in the cal- 
culation, ut\ n is still increasing [a>i n = 148 K (x n = 1.00), 
w ln = 154 K (x n = 0.85), w ln = 155 K (x n = 0.75)]. 

This unambiguously demonstrates that the influence 
of carbon on the low-energy modes, which are crucial for 
the observed strong electron-phonon coupling in MgCNi 3 
(see Ref. IT3) can not be neglected. A similar result was 
already derived from structural investigationsii In order 
to stress this point, A p h was calculated from Eq. 
using S w 4.8 from Ref. 0. Fig. shows a comparison 
of Aph determined by Eq. (JJJ and Eq. J2J, respectively. 
From the phonon density of states a linear dependence 
from the effective carbon content x c is derived (open sym- 
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Table I: Experimental quantities and electron-boson coupling 
constants as derived from the electronic specific heat. 
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Figure 6: Electron-phonon coupling constant vs. effective 
carbon content. Filled symbols: determined from Eq. Q. 
Open symbols: calculated from Eq. (0. The lines are guides 
to the eye. 



the electron-phonon interaction function a 2 (to) , entering 
Eq. ©. 

The present results clearly point to the importance 
of carbon for stabilizing the low-energy Ni-dominated 
phonon modes which are expected to lead to the high 
electron-phonon coupling in this compound^ A sim- 
ilar scenario was established by iJohannes and Pickettl 
who attributed the phonon hardening in carbon-deficient 
MgCNi 3 and the related carbon-deficient ZnCNi3 to 
the low-energy Ni "breathing" mode^fi It seems un- 
likely that the observed strong decrease of A p h is ex- 
plained by the decreasing lattice parameter only, since 
pressure experiments show an increase of T c with in- 
creasing pressure Considering a possible decrease 
of the ED OS, 5 one would expect slightly higher values 
of Aph- However, considering the well known McMil- 
lan formula and the low T c < 2.5 K of the samples 
with x n < 1.00 this effect seems to be negligible. 6 The 
present results may also help to understand recent in- 
vestigations of the carbon isotope effect performed by 
iKlimczuk and Caval^ which is much stronger than ex- 
pected from calculationsii^iii Further experimental and 
theoretical studies of the evolution of the EDOS as well 
as the PDOS of carbon-deficient and isotope-pure sam- 
ples and the behavior of the electron-phonon interaction 
function a 2 (ui) are highly desirable. 
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